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Abstract: We report on the results of numerical simulations of 1 + 1 dimensional SU(N) 
Yang-Mills theory with maximal supersymmetry at finite temperature and compactified on a 
circle. For large N this system is thought to provide a dual description of the decoupling limit 
of N coincident Dl-branes on a circle. It has been proposed that at large N there is a phase 
transition at strong coupling related to the Gregory-Laflamme (GL) phase transition in the 
holographic gravity dual. In a high temperature limit there was argued to be a deconfinement 
transition associated to the spatial Polyakov loop, and it has been proposed that this is the 
continuation of the strong coupling GL transition. Investigating the theory on the lattice for 
SU (3) and SU (4) and studying the time and space Polyakov loops we find evidence supporting 
this. In particular at strong coupling we see the transition has the parametric dependence 
on coupling predicted by gravity. We estimate the GL phase transition temperature from the 
lattice data which, interestingly, is not yet known directly in the gravity dual. Fine tuning in 
the lattice theory is avoided by the use of a lattice action with exact supersymmetry. 
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1. Introduction 

In recent years a series of numerical studies have been undertaken to explore and test conjec- 
tured holographic dualities between supersymmetric gauge theories and supergravity theories 
[|]]. So far these studies have been confined to the case when the super Yang- Mills theory 
is one dimensional and the dual gravitational theory describes the low energy dynamics of 
DO-branes |, @, |, § § |7|, g, § or the M = 4 theory compactified on S 3 x R [[R], [TT], |12| 



In this paper we extend these calculations to the case of N coincident Dl-branes wrapped 
on a spatial circle, which in the decoupling limit is described by a two dimensional maximally 
supersymmetric Yang-Mills (SYM) theory on a circle @, (See @ H, |l|, for some 



recent developments on the lattice along this direction.) This two dimensional Yang-Mills 
system possesses a richer structure at large TV" than its one dimensional counterpart, since 
when one compactifies the spatial direction on a circle, one can construct a new dimensionless 
coupling that can be varied in addition to the temperature. Arguments from a high tempera- 
ture limit and also from strong coupling, using a dual supergravity description, indicate that 
the system should possess an interesting phase structure in the 2d parameter space spanned 
by the temperature and this new coupling in the large N limit |l3|, [18]. A large N transition 
between confined and deconfined phases with respect to the spatial Polyakov line is expected 
which interpolates between the high temperature region and the strongly coupled region. In 
particular for the strongly coupled region the dual Dl-brane system can be described by 
certain black holes in supergravity, with a compact spatial circle. Then arguments from the 
dual gravity model indicate a first order Gregory-Laflamme (GL) phase transition between 
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the black hole solutions localized on the circle and uniform black hole solutions which wrap 



the circle gg g |2l], |22], |TJ, |§ ||, |||. Translating back to the SYM, the dual gravity model 
predicts the parametric dependence of the transition temperature against dimensionless circle 
coupling - a dependence which seemingly cannot be deduced by simple SYM considerations. 
Interestingly, since the relevant gravity solutions have not been constructed yet (analog solu- 



tions are known, but not in the correct dimension [25, 26, pTj), the precise coefficient in this 



relation is not known, and determining it in SYM yields a prediction for the phase transition 
temperature that could be tested in the future when the gravity solutions are constructed - 
a classical but nonetheless rather non-trivial gravitational problem. 

The purpose of this paper is to use Monte Carlo simulation of a lattice formulation of the 
two dimensional SYM theory to investigate its phase structure, focussing on possible large 
N transitions between spatially confined and deconfined phases of the model as revealed by 
behavior of the spatial Polyakov line. In the next section we review the theoretical background 
to the conjectured 2d Yang-Mills/Dl-brane duality when the theories are compactified on a 
circle and describe the expected phase structure in certain limits. 

The usual problems associated with the study of supersymmetric lattice theories are 
avoided by use of new formulations which possess exact supersymmetric invar iances at non- 
zero lattice spacing - the relevant lattice construction is described in section ^. We then 
describe our numerical results, which appear to confirm the expected deconfinement phase 
transition. Furthermore, at strong coupling the position of the observed critical line agrees 
with the parametric dependence on couplings predicted by the dual gravity analysis. In 
particular we estimate the coefficient in this relation and hence derive a prediction for the GL 
phase transition temperature for the dual black holes theory. We believe that this is the first 
time holography has made a detailed prediction for properties of some currently unknown 
non-trivial solutions in classical gravity based on calculations in strongly coupled Yang-Mills. 



2. Theoretical background 

We are interested in studying large N thermal two dimensional maximally supersymmetric (16 
supercharge) SU(N) Yang- Mills theory, in the 't Hooft limit, with coupling A = iV<?y M , with 
the spatial direction compactified. Continuing the theory to Euclidean time, this implies the 
Yang-Mills theory is defined on a rectangular 2-torus, with time cycle size /3, and space cycle 
size R. The fermion boundary conditions distinguish the two cycles, being anti-periodic on 
the time cycle so that f3 has the interpretation of inverse temperature, and periodic boundary 
conditions on the space cycle. The action may then be written as, 

S = jf drdxTx I \l% + \ YlM 1 ' D ^ ~ I ' ^ + fermions ) (2- 1 ) 

Jt * \ i i,J J 

where I, J = 1, . . . , 8 and (j) 1 are the 8 adjoint scalars, and r is the coordinate on the time 
circle, and x the coordinate on the space circle. Since A, f3 and R are dimensionful, it is 
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convenient to work with the two dimensionless couplings, r T = y/~\/3 and r x = \/XR which 
give the dimensionless radii of the time and space circles respectively, measured in units of 
the 't Hooft coupling. We will be interested in the expectation values of the trace of the 
Polyakov loops on the time and space circles, 



N\ 



Tr (Pexp(i<bA T )) ) , P x = 1( Tr (Pexp(i j> A x )) 



(2.2) 



as at large N, these give order parameters for confinement /deconfinement (or center symmetry 
breaking) phase transitions which we will discuss below. 

As discussed in [13, ^] there are several interesting limits for the theory. In the large 



torus limit, 1 <C r x , r T the string theory dual may be described by supergravity. For the weak 
coupling limit, r x ,r T <C 1, or asymmetric torus limits r T <C r x and r x <C r^, we will find the 
dynamics is captured by a lower dimensional YM theory. Let us now review these cases and 
their predictions. 

2.1 Large torus limits and IIB and IIA supergravity duals 

When the torus becomes large in units of the 't Hooft coupling one finds that in certain 
regimes the dual Dl-branes in string theory can be well described by supergravities [|lj as we 
shall now briefly review. Having a supergravity description of the full string theory dual allows 
certain behaviours of the theory to be studied using simple semi-classical gravity reasoning 
which allows powerful predictions to be inferred for the dual SYM. 

The dual IIB string theory is given by the 'decoupling limit' of N coincident Dl-branes 
[|J. This decoupling limit is where one considers finite energy excitations of the Dl-branes 
while taking the limit, gyM = 2^c7 = fi xeci > a ' ~^ 0> where g s is the string coupling and 
a' determines the string tension. Since our Euclidean SYM is defined on a torus, the string 
dual is too, being at finite temperature and having one spatial direction compactified into a 
circle radius R with periodic fermion boundary conditions. 

One finds that for 1 <C r T <C r x this string theory can be described effectively by its 
supergravity sector. String oscillator and winding mode corrections to this supergravity 
description are small in this limit. The IIB supergravity solution describing the thermal 
vacuum is a black hole, carrying electric Dl-brane charge. The Dl-brane charge is string 
like (ie. its field strength tensor is a 3- form), and the appropriate configuration is to take the 
charge to wrap over the compact space circle. The solution preserves translational invariance 
around the space circle direction and is thought to be stable to small perturbations. 

However there is a second supergravity description of the theory which is valid in a partly 
overlapping and partly complementary range 1 <C r T and r x ^ 3 <C r T , obtained by performing 
a T-duality transformation on the compact spatial circle of the IIB string theory [13, 18]. 



Roughly speaking, such a T-duality exchanges winding and momentum modes of the string 
on this spatial circle, and exchanges the IIB string theory for a IIA string theory. In our case 
the N Dl-branes now get exchanged with N DO-branes in the IIA theory. Since the DO-branes 
are point like, rather than string like, they have freedom to distribute their electric charge 
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Figure 1: Figure indicating the regions of coupling space where at large N the dual string 
theory may be approximated by (red) IIB supergravity and (blue) IIA supergravity. In these 
regions, the SYM thermodynamics is dual to the thermodynamics of certain black holes in the 
corresponding supergravity. The IIA region predicts a large N first order phase transition (the 
Gregory-Laflamme phase transition) between black holes localized on the spatial circle, and 
wrapping over the circle. The phase transition is known to occur along the curve = c cr ^r T 
where c cr u is a constant, not yet determined, but known to be order one and c cr n > 2.29. 
The SYM transition is thought to be a deconfinement transition of the spatial Polyakov loop. 

over the circle in various ways - it may be uniformly distributed, non-uniformly distributed, 
or fully localized on the circle, the latter two choices breaking the translational symmetry 
along the space circle direction. It is then a dynamical question which case is preferred. 

It is thought ]l3| that there are 3 types of black hole solution which indeed realize these 3 
choices. The uniform black hole solution exists for all temperatures, but it is known to have a 
dynamical perturbative instability of the Gregory-Laflamme type pq , p9[ for low temperatures 
r\ < 2.29 r T [O]. For higher temperatures it is thought to be dynamically stable. However, 
at a higher temperature than the instability point, so that = c CT n r T for some constant c cr n 
with c cr it > 2.29, the uniform black hole is thought to become globally thermodynamically 
less favored than the localized black hole solution. The actual transition temperature which 
governs the constant c cr a is not yet known, as the localized black hole solutions have not 
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yet been constructed in the correct context to be embedded in the supergravity dual. The 
line = c cr n r T represents a first order thermal phase transition between the uniform and 
localized solutions, with uniform favored for higher temperature r x > c cr u r T and localized 
favored for lower temperature r x < c cr ^r T . We term this the GL phase transition and 
emphasize that this is distinct from the GL dynamical instability. Whilst there is a non- 
uniform black hole solution it is never thermally dominant. For reviews on the GL dynamical 
instability, phase transition and uniform, non-uniform and localized black hole solutions see 

© 0,0- 

According to the duality hypothesis, a Polyakov loop about the time/space circle in 
the Euclidean SYM is computed in the leading large N limit by considering whether a two 
dimensional minimal area surface (the classical string worldsheet) that asymptotically wraps 
the time/space circle exists. If the time/space circle is contractible in the interior of the 
gravity solution, a minimal area solution for the string worldsheet will exist and then the 
correspondence states that P T / X ~ 0(1). However, if the circle is not contractible, there 
cannot exist a minimal surface that gives a finite action for the string worldsheet, and the 
correspondence states that P T / X ~ 0(1/N) and hence P T / X = in the large N limit. It is a 
standard result of Euclidean gravity that black hole solutions have contractible time circles 
in the interior of the solution, and in fact the time circle contracts precisely at the horizon. 
The contractability of the spatial circle however depends on the type of black hole. In the IIB 
supergravity solution the space circle is non-contractible. The IIA uniform (and non-uniform) 
solutions have non-contractible space circles, whereas the localized solution has a contractible 
circle. In fact the eigenvalues of the SYM spatial Polyakov loop (which are phases, and hence 
live on a circle) are thought to correspond to the positions of these DO-branes on the space 
circle in the IIA dual. Hence the GL phase transition can physically be thought of as a thermal 
instability associated with the clumping of DO-branes, breaking the U(l) circle translation 
symmetry. In the large N SYM this symmetry breaking is the spontaneous breaking of center 
symmetry Zn, where for large N, U(l) ~ Z^. 

Let us summarize our predictions for the large torus. We learn that in the IIB regime, 
1 <C r T <C r x , we expect P T ^ but P x = 0. In the IIA regime, where lCr T and r x ^ -C r T , 
we have P T ^ 0, and P x ^ for r x < c cr itr T and P x = for r x > c cr itr T , with c cr u an 
order one constant with c cr a > 2.29. We note that in the regime where both IIA and IIB 
apply, they give consistent results. Thus in the large torus, supergravity regimes, the SYM is 
always deconfined in the time direction, and there is a first order deconfinement /confinement 
transition in the space direction at r x = c cr itr T . 



1 The region r 2 x < or, for 2.29 < a < c cr i t is the region where the localized solution dominates the canonical 
ensemble, but the uniform phase could in principle be constructed as a metastable supercooled state (although 
here we will only be concerned with equilibrium thermodynamics). 



- 5 - 



2.2 Dimensional reduction 

Consider the toy model scalar theory defined on the 2-torus, 

S= \ f drdx {(d^) 2 + 4 ) (2.3) 

A J T 2 

First we change to angular coordinates 9 T = t//3 and 6 X = x/R with unit radius, so 9 T ^ X ~ 
®t,x + 1) and then define the dimensionless scalar variable cf> = {(3R/X) 1 ^ (f>. The action can 
now be written as, 

S = r dO T d6 x U A + ^(dej) 2 + MideJ) 2 ) (2.4) 
Jo V V t V x J 

and we see that the dimensionless couplings r x /rf and v T jr\ determine the masses of the 
non-constant modes of the field <f> on the torus. There are three interesting limits. When 
r x ~ r T <C 1, then the non-constant modes of the scalar become very massive and hence weakly 
coupled and one may integrate these out to arrive simply at the quartic integral governing 
the constant modes. If only 1 <C r x jr\ then the non-constant modes on the time circle are 
weakly coupled and one may integrate these out to obtain the dimensional reduction which 
now lives only on the space circle. Likewise, if 1 r T /r x , one may dimensionally reduce to 
obtain a theory only on the time circle. 

The structure of this toy example is such that precisely the same phenomena occurs with 



the full SYM on a 2-torus, as discussed in |18]. One difference is that due to the anti-periodic 
boundary conditions on the time circle, the Fourier decomposition of the fermions contain 
only non-constant modes in the time direction. Another difference is that under a reduction, 
the constant component of the gauge field in the direction of reduction yields a scalar field, 
similar to the scalars <j> , in the reduced theory. This scalar in the reduced theory corresponds 
to the Polyakov loop about the cycle that has been reduced on. Since the expectation value of 
the eigenvalue distribution of the scalar in the reduced theory will have a non trivial profile, 
this implies that center symmetry is broken in the Polyakov loop about the reduced cycle 2 . 

There are again 3 regimes. For r x ~ r T <C 1 one may reduce on both time and space to 
just give the zero modes of the theory, and arrive at a bosonic Yang-Mills matrix integral, 
since in reducing on the time circle one loses the fermions which have no zero modes. Such a 
reduction indicates that in this limit, the two dimensional SYM should have P T ,P X ^ 0. 

For r\ <C r T the theory may be dimensionally reduced on the space circle to give the 
thermal super symmetric matrix quantum mechanics living on the time circle with radius (5. 
The spatial Polyakov loop is then given in terms of one of the 9 scalars of the BFSS model, and 
since these scalars have localized eigenvalues, the two dimensional SYM should be deconfined 
in the space direction with P x ^ 0. This theory is precisely the BFSS theory |34j|, and 
recently this has been numerically simulated in the 't Hooft limit ||, f|, 0, and indeed, the 



2 This is to be contrasted with the Eguchi-Kawai reduction |33| ] where quite the opposite occurs; one can 
only reduce on a direction if center symmetry is unbroken. 
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Figure 2: Figure showing the regions of coupling space where the SYM may be dimensionally 
reduced on the time and/or space circles. The blue region indicates the region where reduction 
on the space circle gives a good approximation, yielding a supersymmetric quantum mechanics 
theory, the BFSS model. The red region indicates where reduction on the space circle to a 
bosonic quantum mechanics (BQM) is a good approximation. This latter reduction predicts 
a large ./V deconfinement phase transition in the spatial Polyakov loop for r 3 = 1.35r T and 
this curve is shown. 



results obtained are consistent with the theory always being deconfined, so P T 7^ 0. The 
coupling of this quantum mechanics is given by r T /( r a:) and when this is large we know 
from our arguments above that we are in a regime where a dual IIA supergravity description 
exists, and the dynamics is given by the localized black hole solution which is indeed consistent 
with P T ,P X ^ 0. 

For r 3 <C r x one may again perform a dimensional reduction, now on the time circle. Thus 
in the two dimensional theory we expect P T 7^ 0. Since there are no fermion zero modes on 
the time circle, the resulting one dimensional theory is a bosonic quantum mechanics (BQM) 
defined on a circle radius R and with dimensionless coupling r 3 /r T . Numerical Jl3| , |l8| , |35| ] 
and analytic study [36] indicates that this theory has a large N confinement /deconfinement 
transition at r 3 /r T ~ 1.35 of second order. There is also thought to be a third order Gross- 
Witten [37, ^] transition very nearby at r^./r T ~ 1.49 [35, 36]. 
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2.3 Expectations for large N phase diagram 

We conclude by putting together the above discussions. The simplest picture is then that the 
Gregory-Laflamme first order phase transition, r x = c cr nr T for 1 <C r T (recall c cr n > 2.29), 
and the second order transition <C r x and r x = 1.35 r T in the time reduced BQM are two 
ends of the same spatial Polyakov loop confinement /deconfinement phase transition line. At 
some point in-between the order presumably changes, and here the new third order Gross- 
Witten phase transition emerges although this is not measured by centre symmetry breaking, 
but by more detailed information about the spatial Polyakov loop eigenvalue distribution. It is 
interesting [ 



13 , 23] that the new phase at small r x also exists for 1 <C r x in the form of non- 



uniform IIA black strings, but unlike at weak coupling, these are never thermally dominant 
in the IIA supergravity region. In figure |3| we summarize the expected phase diagram for the 
spatial confinement/deconfinement transition. 




t r T , c crit >2.29 



P_t + everywhere 



Figure 3: Cartoon of the expected large N, spatial Polyakov loop deconfinement transition line 
in coupling space. Pictured is the simplest possibility, namely that the spatial deconfinement 
transition interpolates between the strong coupling Gregory-Laflamme transition parametric 
behavior ~ r T , and the high temperature reduction deconfinement transition behavior 
rt = 1.35r T . 



3. Supersymmetric lattices for super Yang-Mills 

In recent years there has been significant progress in the formulation of lattice theories which 
possess exact supersymmetry - see |3^, £l0| and references therein. The chief advantage of these 
approaches is that they significantly reduce, or in many cases eliminate, the need to fine tune 
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the couplings in the lattice theories to approach the target continuum super symmetric field 
theory as the lattice spacing is sent to zero. The case of M = 4 super Yang-Mills is particularly 
interesting and the corresponding lattice theory was first constructed using orbifold methods 



in [41|. Subsequently, it was realized that the same theory could be obtained using a carefully 
chosen discretization of a topologically twisted version of the continuum theory ^2|. [43| . 

The continuum twist of M = 4 that is the starting point of the twisted lattice construction 
was first written down by Marcus in 1995 H] although it now plays an important role in the 



Geometric-Langlands program and is hence sometimes called the GL- twist [45|. This four 
dimensional twisted theory is most compactly expressed as the dimensional reduction of a 
five dimensional theory in which the ten (one four component gauge field and six scalars) 
bosonic fields are realized as the components of a complexified five dimensional gauge field 
A a , a = 1 ... 5 while the 16 single component twisted fermionic degrees of freedom are realized 
as the 16 components of a Kahler-Dirac field (j],tp a ,Xab) The appearance of a scalar 

fermion after twisting is crucial - it implies the existence of a nilpotent supersymmetry which 
will be preserved in the lattice theory. Its action on the continuum fields 3 is given by 

Qi> a = o 
QA = o 

Q Xab = -F a b 
QV = d 

Qd = 

The site field d is a non-dynamical field that is included to close the Q-algebra and is subse- 
quently integrated out of the final lattice action. Furthermore, the action of this theory can 
be written as the sum of two terms; a Q-exact piece of the form 

S = J-S / Tr {xabTab + Vpa,T> a ] - Ivd) (3.1) 



9ym 

and an additional Q-closed term 



2 



If — 

•^closed = ~ 2 / Tr CrnnpqrXqr'DpXrnn (^-2) 

^9ym J 

The super symmetric invariance of this term then relies on the Bianchi identity 

^-mnpqr^p-F qr — (^ - ^) 

Discretization of this theory proceeds straightforwardly; (Complex) continuum gauge fields 
are represented as complexified Wilson gauge links Unix) = e" 4 ' 1 ^ living on links e-, jjb = 



3 We assume an anti-hermitian basis for all fields which take their values in the adjoint representation of 
the SU(N) gauge group. 



- 9 - 



1 ... 4 of a hypercubic lattice. The field U§ is an exception to this and is placed on the body 
diagonal of the hypercube corresponding to a relative position vector e$ = (—1, —1, —1, —1). 
Notice that Yla=i e a = which is important to show gauge invariance of the lattice action. 
These fields transform in the usual way under U (N) lattice gauge transformations eg: 

U a (x) -> G(x)U a {x)G\x + e a ) (3.4) 

Supersymmetric invariance then implies that ip a { x ) nve ° n the corresponding links and trans- 
form identically to lA a {x). The scalar fermion n(x) is clearly most naturally associated with 
a site and transforms accordingly 

rj(x) G{x)n{x)G\x) (3.5) 

The field \ab is slightly more difficult. Naturally as a 2- form it should be associated with a 
plaquette. In practice we introduce diagonal links running through the center of the plaquette 
and corresponding to the vector e ab = e a + e b and choose Xab{x) to lie with opposite orienta- 
tion along those diagonal links. This choice of orientation will again be necessary to ensure 
gauge invariance. The scalar lattice supersymmetry transformation is identical to that in the 
continuum after the replacement A a — > U a - Most importantly it remains nilpotent which 
means that we can guarantee invariance of the Q-exact part of the lattice action by replacing 
the continuum fields by their lattice counterparts. Of course to do so necessarily requires a 
prescription for replacing continuum derivative operators by gauge covariant finite difference 
operators. The following expressions are used: 

vMf b = U a (x)f b (x + ea) - h(x)U a (x + e b ) (3.6) 
V { ~ ] f a = fa{x)U a {x) -U a {x- e a ) f a (x - e a ) (3.7) 

Notice that this definition reduces to the usual adjoint covariant derivative in the naive con- 
tinuum limit and furthermore guarantees that the resultant expressions transform covariantly 
under lattice gauge transformation. The lattice field strength is then given by the gauged 
forward difference J- ab = T>^U b and is automatically antisymmetric in its indices. Further- 
more it transforms like a lattice 2-form or plaquette and hence yields a gauge invariant loop 
on the lattice when contracted with the plaquette fermion Xab- Similarly the covariant dis- 
crete divergence appearing in T>^ a ^hi a transforms as a 0-form or site field and hence can be 
contracted with the site field rj to yield a gauge invariant expression. 

This use of forward and backward difference operators guarantees that the solutions of the 
theory map one-to-one with the solutions of the continuum theory and hence fermion doubling 
problems are evaded [|46[| . Indeed, by introducing a lattice with half the lattice spacing one 
can map this Kahler-Dirac fermion action into the action for staggered fermions [47]. While 



4 A better choice in four dimensions is the A\ lattice which retains a higher point group symmetry than the 
hypercubic lattice. See [^9| for details. It is not necessary for two dimensions and indeed would complicate 
the calculation of Polyakov lines. 
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the super symmetric invariance of the Q-exact term is manifest in the lattice theory it is not 
clear how to discretize the continuum Q-closed term. Remarkably, it is possible to discretize 
( |3.2[) in such a way that it is indeed exactly invariant under the twisted supersymmetry 

^closed = 2 ^ ^ Tr CmnpqrXqr ( x "t~ &m ~\~ &n ~\~ ^p)^p Xmn( x ~l~ (3-8) 

°9ym x 

and can be seen to be super symmetric since the lattice field strength satisfies an exact Bianchi 
identity @. 

f-mnpqrDp -F qr — (3-9) 

Putting all these elements together we arrive at the supersymmetric lattice action |39| 
S = S closed + -L Tr Ulb^ab + \ (v { - ] U a ) 2 - XabV { + ] ^ b] - rfD { ~ Va) (3-10) 

where we have taken the Q-variation and integrated out the auxiliary field d. To reiterate; 
this action is gauge invariant, free of doublers and possesses the one exact supersymmetry 



given in (3.1). 

Finally to obtain a two dimensional theory we perform a simple dimensional reduction 
along two lattice directions using periodic boundary conditions. The resultant lattice action 
corresponds in the naive continuum limit to the target Q = 16 YM theory in two dimensions. 
In this limit its exact supersymmetry is enhanced to correspond to 4 continuum supercharges 
corresponding to the four scalar fermions that now appear in the dimensionally reduced theory 

m 

We will be interested in this theory large N limit with 't Hooft coupling A = Ngy^.j- The 
lattice theory is then governed by the coupling k = where L and T denote the number 
of lattice sites in the spatial and temporal directions. 

We have used periodic boundary conditions for the fields on the remaining spatial circle 
and anti-periodic boundary conditions for fermions in the temporal direction in order to 
access the thermal theory. Simulations were carried out using the RHMC algorithm which 
is described in detail in [49[ |. It has been shown that the existence of a noncompact moduli 
space in the theory renders the thermal partition function divergent Q. In order to regulate 
this divergence we have additionally introduced a mass term for the scalar fields appearing 
in the lattice action with a dimensionless mass parameter m = m p h ys /?- 



S 



9 

m \ ^ 



m 2 / j 

9ym t 



(3.11) 



The form of this term is effective at suppressing arbitrarily large fluctuations of the exponen- 
tiated scalar fields and reduces to a simple mass term for small fluctuations characterizing 
the continuum limit. Notice that this infrared regulator term breaks supersymmetry softly 
and lifts the quantum moduli space of the theory. We have performed our simulations for a 
range of the parameter m in order to allow for an extrapolation m — > 0. 
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4. Simulation results 
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Figure 4: Spatial and temporal Polyakov lines (P x and P T ) against dimensionless time circle 
radius r T for maximally supersymmetric SU(3) Yang-Mills on 2 x 8 and 3 x 12 lattices using 
different values of the infrared regulator m. 



In this section we present our numerical results. We have focused on the Polyakov lines 



for both the thermal and spatial circle. These are defined in the usual way 



P x = -( 
x N\ 



Tr Ii L a ;},U ax ), P T = 1/ Tr ) , (4.1 



N\ 



where we have extracted the unitary piece of the complexified link to compute these 
expressions. We have evaluated the spatial and temporal Polyakov lines as a function of r T 
for two different lattices with the same aspect ratio, a 2 x 8 lattice and a 3 x 12 lattice, for 
N = 3 and with values of the infrared regulator m = 0.05,0.10 and 0.20. The use of two 
different lattices with the same aspect ratio allows us to test for and quantify finite lattice 
spacing effects. We have performed simulations for values of the dimensionless time circle 
radius in the range 0.02 < r T < 1.0. The results are shown in figure [|. 

Notice that the temporal Polyakov remains close to unity over a wide range of r T . This 
indicates the theory is (temporally) deconfined and is consistent with expectations for the 
limits discussed in section ^ - the asymmetric torus limits, and the strong coupling regions 
where there is a dual supergravity description in terms of black holes. However, the spatial 
Polyakov line has a different behavior taking values close to unity for small r T while falling 
rapidly to plateau at much smaller values for large r T . It is tempting to see the rather rapid 
crossover around r T ~ 0.2 as a signal for a would be thermal phase transition as the number 
of colors is increased. This conjecture is seen to be consistent with the data; in figure |5] we 
show the Polyakov lines for JV = 2,3,4 on 2 x 8 lattices as a function of r T . The plateau 
evident at large r T falls with increasing N and the crossover sharpens. This is consistent with 
the system developing a sharp phase transition in the large N limit. 

In the data shown here we see very little dependence of our results on the scalar mass. 
Indeed for the length of Monte Carlo we were able to perform it appears that m can be set 
to zero for r T < 2 without fear of encountering the thermal divergence discussed in [|| . This 
stability in the scalar sector can be seen in figure || which shows the Monte Carlo time series 
for the eigenvalues of U^U^ ~ e 2 ^ at two different r T 's with dimensionless mass parameter 
m = 0.05 and gauge group SU(3). There is no evidence of a divergence over thousands of 
Monte Carlo sweeps. Furthermore, one sees that the eigenvalues of the scalar fields (rendered 
dimensionless using the lattice spacing) cluster with small separation for this range of r r . 5 

We have however observed that the m = model does exhibit the same thermal instability 
observed in the case of supersymmetric quantum mechanics for sufficiently low temperature 
r T >> 1 in agreement with the general arguments given in 

Putting together several lattice aspect ratios for N = 3, 4 we can plot the spatial Polyakov 
loop as a function of r s and r T where data is available. This is done in figure fj]. The three 
contours P x = 0.4,0.5,0.6 are shown. We see that the contours for SU(4) are closer together 
than those for SU (3), as we expect for a large N transition. From this data we can try to assess 
where the large N transition in P x may occur. In the detailed studies of the dimensionally 



reduced bosonic quantum mechanics 13], it was found that the large N transition occurred 



5 We thank Masanori Hanada for pointing out the significance of this result which differs from the situation 
reported in Mm 
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Figure 5: Plot of the absolute values of the spatial and temporal Polyakov lines {P x and P T ) against 
the dimensionless time circle radius r T for maximally supersymmetric SU(N) Yang-Mills on a 2 x 8 
lattice for N — 2,3, 4, using the value of the infrared regulator m = 0.10. 



very close to P x ~ 0.5. Thus from the contours of the SU (3) and SU (4) data we could take the 
P x = 0.5 curves to give an estimate for the large N phase transition line. Another estimate is 
to plot the function /„ = P x (SU(n)) — P x (SU(n — 1)) which measures the difference between 
the Polyakov lines for SU(n) and SU(n— 1). At strong coupling where we expect the large N 
transition is first order, the simplest situation is to have f n < in the confined region (where 
P x = for TV — > oo), and correspondingly f n > in the deconfined region as n — > oo. Then 
plotting the boundary of the positive (or negative) region of calculated from our data also 
gives an estimate of the critical line. Neither method can give a precise determination, and 
they should not be considered as a replacement for calculations at larger than we have been 
able to reach here. However, in the absence of such large N data we plot the P x = 0.5 contours 
for 5f7(4) and SU(3) in figure 0, and in addition the region where is positive. We note that 
the SU(3) and SU(4:) P x = 0.5 contours are remarkably consistent with each other, which 
provides evidence that they are indeed a reasonable approximation to the large N transition 
curve. Whilst the data is rather noisy, and hence the positive region has 'holes' in it, 
the function is positive only to left of the P x = 0.5 curves, and furthermore, extends right 
up to these curves. The curve = 3.5r r is plotted on this graph and matches the contours 
P x = 0.5 and the boundary of the positive region very well in the strong coupling region. 
We take this to indicate that the gravity prediction for the parametric behavior = c cr ^r r 
is consistent with our data, and we have estimated c cr u — 3.5 which indeed obeys the gravity 
prediction that c cr u is order one and c cr u > 2.29. Furthermore we see that the contours 
P x = 0.5 also appear to be consistent with the high temperature prediction = 1.35r T as 
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Figure 6: Plots of the average scalar eigenvalues against Monte Carlo configuration time step, 
for iV = 3ona2x8 lattice with r T = 0.5 and 1.0. Note that the spread between eigenvalues 
reduces as r T is decreased. We have used the dimensionless mass parameter m = 0.05. 



well. 

The value of the ratio a = c cr #/2.29 gives the ratio of the GL thermal phase transition 
temperature to the GL dynamical instability temperature (the minimum temperature to 
which uniform strings can be supercooled), so a = Tgl phase /TGLinstab- Whilst the GL 
instability temperature is known |l| (corresponding to the behavior r\ = 2.29r T at strong 
coupling), the GL phase transition temperature is not known in the gravity theory as the 
localized solutions have not been constructed. In fact the near extremal DO-charged black 
holes are simply related to vacuum solutions of pure gravity with R 1,8 x S 1 asymptotics 
[ p"ii| , |23[ . Such localized black hole solutions have been constructed for asymptotics 1R 1,3 x S 1 
and IR 1,4 x S , using numerical techniques []50], 25, p7[ ]. Extending these methods to the case of 
interest here, R 1 ' 8 x S , is obviously an interesting future direction. It is worth emphasizing 
that whilst finding localized solutions in the gravity theory only involves solving the classical 
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Figure 7: Plot of contours of the expectation of the spatial Polyakov line P x over the r x ,r T plane. 
The left frame shows SU(3), and the right 5[/(4). The three contours plotted are 0.4, 0.5, 0.6, and the 
simulation data collates and interpolates runs made on lattices 2 x 16, 2 x 8, 3 x 8, 4 x 4 and 4x8 
therefore giving a variety of aspect ratios r T /r x . 



Einstein equations, in practice even phrasing the Einstein equations in a manner amenable to 
numerical solution has presented a challenge [27] and then solving the resulting coupled partial 
differential equations is a serious numerical undertaking. 6 Our lattice estimation a ~ 1.5 
provides a prediction for the thermal behavior of the gravity solutions. To our knowledge, 
this is the first time a prediction about the properties of non-trivial classical gravity solutions 
has been made from the Yang-Mills side of a holographic correspondence. 



5. Conclusions 



We have used a supersymmetric lattice action to study the strongly coupled dynamics of two 
dimensional maximally supersymmetric SU (N) Yang-Mills theory at finite temperature and 
compactified on a circle for a range of N = 2, 3, 4. In particular we have focused on the 
spatial Polyakov line as an order parameter for a large N deconfining phase transition. Our 
simulations are consistent with the existence of a single transition curve in the 2d parameter 
space spanned by the two dimensionless couplings, r x ,r T which give the size of the thermal 
and spatial circle in units of the YM coupling. 

At high temperature, <C r x , the simulations are consistent with the previously pre- 
dicted behavior that the transition curve goes as = 1.35r T . At strong coupling, 1 <C r T , the 

6 Such solutions can be constructed perturbatively J5l], |5^, Q in a small radius limit (compared to the 
circle size) but the GL phase transition occurs for black holes with radius of order the circle size, and hence it 
is unclear how accurate perturbative methods are for a prediction of Tgl phase/TGL instab- 
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Figure 8: Plot showing a superposition of the P x = 0.5 contours for SU(3) and SU(4) as dashed 
black lines. Also shown is the region (blue) where the SU(4) loop P x is greater than the SU(3) 
loop, which is expected to estimate the large N deconfined region for a first order transition (which 
gravity suggests at strong coupling). 'Holes' in this blue region are due to statistical errors. We see the 
boundary of this region (ignoring 'holes') matches well the P x = 0.5 contours, and represents our guess 
for where the large N transition resides. This figure should be compared to the previous figure || giving 
a sketch of the expected phase structure. Plotted on the figure is the high temperature prediction for 
the transition (r x = 1.35r T , red curve). We note that the estimated large N transition curve fits well 
both this high temperature prediction and also the strong coupling dual gravity predicted parametric 
behavior r 2 z = c cr itr T . Our data suggests c cr a — 3.5 (plotted as blue curve), which obeys the constraint 
from gravity c cr u > 2.29. 

transition is conjectured to be the holographic dual of a first order Gregory-Laflamme phase 
transition, with the transition curve going a r%, = c cr itr T , with c cr a an order one constant 
obeying the constraint c CT n > 2.29. Our simulations are consistent with this parametric be- 
havior, and we have used the N = 3, 4 data to estimate the position of the large N transition, 
determining c cr u — 3.5. This gives the ratio of the Gregory-Laflamme phase transition and 
dynamical instability temperatures to be Tahphase/TGLinstabiUty — 1-5- Since the dual local- 
ized black hole solutions have not been constructed, this constitutes a prediction for these 
non-trivial gravity solutions, which hopefully will be tested by their construction in the near 
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future. 
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